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Abstract 


We use Monte Carlo simulation to study the phase diagram of three-dimensional dynamical 
triangulations with a boundary. Three phases are indentified and characterized. One of these 
phases is a new, boundary dominated phase; a simple argument is presented to explain its 
existence. First-order transitions are shown to occur along the critical lines separating phases. 

Dynamical triangulations with a boundary term 

Dynamical triangulation models arise from simplicial discretizations of continuous Riemannnian 
manifolds. A manifold is approximated by glueing together a set of equilateral simplices with 
fixed edgelengths. This glueing ensures that each face is shared by exactly two distinct simplices 
- the resultant simplicial lattice is called a triangulation. In the context of Euclidean quantum 
gravity it is natural to consider a weighted sum of all possible triangulations as a candidate for 
a regularized path integral over metrics. Physically distinct metrics correspond to inequivalent 
simplicial triangulations. This prescription has been shown to be very successful in two-dimensions 
(see, for example, 0). 

Most analytic studies and almost all numerical work done so far has been restricted to compact 
manifolds like the sphere. In this paper we develop techniques that allow us to extend numerical 
studies to simplicial manifolds with boundaries. Specifically, we study the 3-disk created by insert¬ 
ing an S 2 boundary into a triangulation of the sphere S 3 . This allows us to compute an object 
which is the simplicial equivalent of the ‘wavefunction of the Universe’ Q: 
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The functional integral over 3-metrics g is restricted to those with 2-metric h on the boundary. 
The simplicial analog is simply 


^(T 2 ) = Y j e~ S ^ ) 

t 3 


( 2 ) 


Thus the probability amplitude for finding a particular 2-triangulation T 2 is obtained by counting 
(with some weight) all 3-triangulations T 3 which contain T 2 as their boundary. A natural lattice 
action Sl can be derived from the continuum action by straightforward techniques [lOj. It contains 
both the usual Regge curvature piece familiar from compact triangulations together with a boundary 
term. The boundary term arises from discretization of the extrinsic curvature of the boundary 
embedded in the bulk. In three-dimensions the curvature is localized on links. If Lm denotes 
the set of links in the bulk of the 3-triangulation (excluding the boundary) and Lqm those in the 
boundary the action can be written 


Sl = «i I ~ an h) + 51 “ an h) I ( 3 ) 
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The quantity a = arccos(l/3) and rih is the number of simplices sharing the link (hinge) h. 
Typically Sl will also contain a bulk cosmological constant that can be used to tune the simulation 
volume. The resultant action can be rewritten in the form 


S b = -k 0 N 0 + k 3 N 3 + K b N b 2 (4) 

where N b is the area of the boundary. Here, k 3 is used to tune the volume of the system. We are 
thus left with a two-dimensional phase space parameterized by kq and K b conjugate to the number 
of vertices and the number of boundary triangles. It is trivial to generalize this to, for example, 
four-dimensions. The partition function for the system is then 

Z = Y,e~ Sb (5) 

T 

where the sum is over triangulations, T. 


Various other extended phase diagrams have been studied for three-dimensional dynamical trian¬ 
gulations including adding spin matter [13, ||, adding gauge matter [12, £|, and adding a measure 
term [^J], Much of this work was motivated by the desire to find a continuous phase transition. 
No such transitions have been found. 


Simulation 


Our simulation algorithm is an extension of the algorithm for compact manifolds in arbitrary 
dimension described by Catterall Q. Consider the environment of any vertex in a D-triangulation 
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- it is composed of simplices making up a trivial D-ball. The boundary of this D-ball is just 
the sphere S( D ~ 1 ' 1 . A boundary with the topology of S can thus be created in the original 
triangulation by removing these simplices. If the original triangulation corresponded to the sphere 
S D the topology of the new triangulation is that of a D-disk. 


In practice we simulate a compact manifold with one ‘special’ vertex. This vertex and all the 
simplices sharing it are ignored during any measurement. In this way every triangulation of our 
marked sphere S D is in one-to-one correspondence with a triangulation of the D-disk. Notice that 
the usual compact manifold moves applied to all simplices (including those sharing the special 
vertex) will in general change the boundary of the D-disk. Indeed these moves are ergodic with 
respect to the boundary. Furthermore, the proof that these moves satisfy a detailed balance relation 
goes through just as for the compact case. The one extra restriction is simple — one must never 
delete the special vertex. With this trick we can trivially extend our compact codes to the situation 
in which a S^ D ~^ boundary has been added. We are merely simulating a compact lattice with an 
action which singles out a special vertex and its neighbour simplices. This contrasts with the set 
of additional boundary moves used by Adi et al Q for simulations in two-dimensions. 


Measurements do not include the special vertex or any simplices connected to it. For example, 
the number of D-simplices in the system with boundary is the number of D-simplices in the whole 
simulation minus the number of D-simplices sharing the special vertex. The size of the boundary 
is simply the number of D-simplices sharing the special vertex. 


We have used the Metropolis Monte Carlo [111 scheme with usual update rule: 

^(accept move) = min{e ~ ASb , 1} (6) 

and in this way we explore the space of unlabeled triangulations with the action Sj, (equation [| for 
three-dimensions). 


Checks in two-dimensions 

In two dimensions we tested our simulation code at small volumes by comparing with hand calcu¬ 
lated amplitudes for small disks. We label disk configurations by the number of triangles and the 
boundary length: (N 2 ,N%). We calculated the ratios of amplitudes for disks (1,1):(2,4):(3,3):(3,5) 
to be 1:1.5:1:3. Our simulation gave 1.01:1.52:1:3.04 from a sample of 1 million disks with volume 
1-3. This test was extended up to volume 5 disks, also showing good agreement. 

We also tested our simulation code by comparing results for two-dimensions of Adi et al [;Ij. All our 
results agree within the statistical errors. Table [I] shows a comparison of the results for a selection 
of lattice sizes. 
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a 2 

(Af) (Adi et al) 

(N±) (this work) 

50 

39.83(5) 

39.88(4) 

100 

78.05(5) 

78.09(4) 

200 

154.54(6) 

154.57(6) 

400 

307.62(8) 

307.53(8) 

800 

613.4(1) 

613.4(1) 

1600 

1225.1(1) 

1225.1(2) 

3200 

2448.7(2) 

2448.8(3) 

6400 

4895.6(3) 

4895.4(13) 


Table 1: Estimates of (iVj’), the expectation value of the boundary size (length), for two-dimensional 
manifolds of various sizes from Adi et al Q, and from this work. 


Phase diagram 


We performed a set of simulations in three-dimensions with action of equation |4]. In all runs K 3 
was used to tune the nominal system volume, A 3 , to 2000 for each given kq and Kb- 

In three-dimensions there are just 4 types of move: vertex insertion, vertex deletion and exchange 
of a link with a face (two moves: link to face or face to link). Where these moves take place 
on sections of the triangulation involving the special vertex we take care to count changes in the 
numbers of simplices inside and outside of the boundary but otherwise the moves are the same as 
for the bulk. Series of runs varying either kq or Kb were made and the vertex susceptibility used to 
search for phase transitions. We define the vertex susceptibility, x, to be normalized with respect 
to the number of 3-simplices: 

X = T«JV 0 2 ) - (No) 2 ) (7) 

The points shown in figure [I] are taken from the positions of peaks in the vertex susceptibility. 

In figure || there are three phases which we characterize as: phase 1 - crumpled, minimal boundary; 
phase 2 - branched-polymer, minimal boundary; and phase 3 - boundary dominated. In phases 1 
and 2 the boundary is simply 4 triangles (2-simplices) connected to form a tetrahedral hole. The 
system is essentially like a compact manifold with one marked 3-simplex — the tetrahedral hole. 
In phase 3 the boundary is large — typically a substantial fraction of the bulk volume. 
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Figure 1: Phase diagram for 3-dimensional dynamical triangulation with a boundary. All points 
have error bars in either Kb or kq, most cannot be seen because they are smaller than the symbols. 
Nominal simulation volume, JV 3 = 2000. 



Figure 2 : Sample of vertex susceptibility data for different values of the boundary coupling 

constant, Kb- 
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Figure 3: Number of vertces, (No), as a function of /to and K b . Nominal simulation volume, 
IV 3 = 2000. Note that we see three distinct areas with different values of (No): the boundary 
dominated phase (small k b , large /to) with large (No), the crumpled phase (small kq ) with small 
(No), and the branched-polymer phase (large n b and kq) with intermediate (No). 


Simple argument for boundary dominated phase 

Here we argue that the boundary dominated phase can be explained by considering an effective 
action written in terms of the boundary size. We show that in certain circumstances a large 
boundary will decrease this action. Otherwise one of the minimal boundary phases will be favored. 

Consider the action: 

S b = -koN 0 + k b N b (8) 

We ignore the volume term as this is kept fixed during the simulation. If we note that the boundary 
is itself a 2 -sphere then we know that: 

N b = 2 (N b - 2) (9) 

and 

No = Nq + IVg (10) 
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where Nq is the number of vertices, Nq is the number of vertices on the boundary, Nq is the number 
of internal vertices, and N% is the number of 2-simplices (triangles) on the boundary. We may thus 
rewrite the action: 

S b = -«o K + ( 2n b - k 0 )N* (11) 

If we now consider Nq fixed and note that the number of manifolds with boundary size N% is 
governed by an exponential factor e K t > N 2 = e 2K b N 0 , where is a new constant, we may then write 
an effective action for the number of boundary vertices: 

S e ff ~ (~«o + 2(«6 - k%))NI (12) 

The presence of small or large boundaries is then determined by the sign of this action. We thus 
expect the phase transition at kq = 2(k& — k%) which is in good agreement with what we see. 


Order of transitions 
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Figure 4: Time series showing the boundary size (IVf) during simulation. The upper plot is at the 
transition between the crumpled and boundary dominated phases (fco = —0.423, k b = 0). The lower 
plot is at the transition between the branched-polymer and boundary dominated phases (fco = 5, 
k b = 2.43). Nominal simulation volume, IV 3 = 2000, and time is in units of IOOIV 3 attempted 
updates. 



0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 

Simulation time, t 



1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 

Simulation time, t 


7 




Simulations of compact manifolds in three and four-dimensions are known to have a first-order 
phase transition between crumpled and branched polymer phases (3d Q, 4d || ||]). Our Monte 
Carlo time series show strong bistability on all three phase boundaries (see figure £|). We take this 
to indicate that all three phase transitions are first-order. 


Concluding remarks 


We have demonstrated an arbitrary dimension algorithm for simulating dynamical triangulations 
with a boundary. This has been tested against known results in two-dimensions and used to map 
the phase diagram in three-dimensions. 

We have identified three phases in three-dimensional dynamical triangulations with a boundary 
and mapped the boundaries within the range of couplings — 1 < kq < 5 and —0.5 < /ct, < 4. 
The observed phases include the crumpled and branched-polymer phases seen in triangulations of 
compact manifolds, and also a new, boundary dominated phase. The existence of this phase, and 
the shape of the phase boundary on the kq-K b phase diagram, is predicted by a simple argument. 
Obvious bistability in the time series at the phase transitions indicates that all transitions within 
the range of couplings studied are first-order. 
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